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Abstract 
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. We review some recent works on the post-Newtonian theory of slowly-moving (post- 

, ^ ' Newtonian) sources, and its application to the problems of dynamics and gravitational 

^ I radiation from compact binary systems. Our current knowledge is 3PN on the center- 

. of-mass energy and 3.5PN on the gravitational-wave flux of inspiralling compact bi- 

' naries. We compute the innermost circular orbit (ICO) of binary black-hole systems 

and find a very good agreement with the result of numerical relativity. We argue that 
the gravitational dynamics of two bodies of comparable masses in general relativity 
does not resemble that of a test particle on a Schwarzschild background. This leads 
us to question the validity of some "Schwarzschild-like" templates for binary inspiral 
which are constructed from post-Newtonian resummation techniques. 
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Recent years have shown a tremendous revival of interest in an Old Lady: the post-Newtonian approxima- 
tion (or expansion when the speed of hght c — > +oo), which is surely the most important technique in the 
arsenal of General Relativity for drawing firm predictions for the outcome of experimental facts related 
to gravitation. The post-Newtonian approximation is ideally suited for describing the adiabatic phase of 
the very interesting astrophysical systems known as inspiralling compact binaries. These systems consti- 
tute our only known-to-exist (for sure) source to hunt for in the current network of laser-interferometric 
detectors of gravitational waves, today composed of the large-scale interferometers VIRGO and LIGO, 
, and the medium-scale ones GEO and TAMA. 

^ ' Two compact (i.e. gravitationally-condensed) objects — neutron stars or black holes — orbit an in- 

ward spiral, with decreasing orbital radius r, decreasing orbital period P, and increasing orbital frequency 
w = The inspiral is driven by the loss of energy associated with the gravitational- wave emission, or, 
equivalently, by the action of radiation forces. There is a long phase of adiabatic inspiral, with associated 
dimensionless adiabatic parameter 

The order of magnitude of the right-hand-side is that of the radiation-reaction force, namely ~ 1/c^ 
or equivalently 2.5PN order beyond the standard Newtonian acceleration (following the usual post- 
Newtonian jargon). The binary's dynamics being essentially "aspherical" , inspiralling compact binaries 
are strong emitters of gravitational radiation. 

The main point about the theoretical description of inspiralling compact binary is that a model made 
of two structureless point-particles, characterized solely by two mass parameters mi and m2 (and possibly 
two spins), is sufficient. Most of the non-gravitational effects usually plaguing the dynamics of binary star 
systems: a magnetic field, an interstellar medium, etc., are dominated by gravitational forces. However, 
the real justification for a model of point particles is that the effects due to the finite size of the compact 
bodies are small. In particular, the tidal interactions between the two compact objects are expected to 
play a little role during most of the inspiral phase; the mass transfer (in the case of neutron stars) does 
not occur until very late, near the final coalescence. Thus the inspiralling compact binaries are very clean 
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systems, essentially dominated by gravitational forces during most of their life. This is the reason why 
these systems are so interesting for relativity theorists! 

Our fascination for the old post-Newtonian lady is that it was recognized that improved waveform 
modelling is crucial for constructing efficient templates for searching and measuring the gravitational 
waves from inspiralling compact binaries in the LIGO/.../TAMA network. By improved we mean in fact 
much improved with respect to what is known from the "Newtonian" approximation, which corresponds 
in the case of the radiation field to the Einstein quadrupole formalism. And one needs much improved 
modelling because the orbital motion of the inspiralling compact binary — which is responsible for the 
gravitational-wave emission — becomes very relativistic when the binary reaches the so-called Innermost 
Circular Orbit or ICO (defined below in this article). At that point the orbital velocity is of the order of 
50% of the speed of light. The price to be paid for applying the post-Newtonian approximation is that we 
must go to very high post-Newtonian order. The status of the field nowadays is the 3PN approximation 
(i.e. — or even better the 3.5PN one — , which is likely to be sufficient for practical purpose (at 

least in the case of neutron stars binaries). 

After the two objects have passed the ICO they will plunge together and merge to form a single black 
hole, which will subsequently settle down into a stationary configuration, by emission of gravitational 
waves in the quasi-normal mode channel. Of course we should not expect the post-Newtonian approx- 
imation (and the description of the compact objects by point-particles) to be valid after the ICO, and 
we shall have to replace this model by a fully relativistic numerical computation of the phase of plunge 
and merger of two black holes T". But, up to the point of the ICO, the post-Newtonian approximation, 
when carried out to 3PN order, is physically valid, and is probably very accurate. This last point may 
appear to be a little surprising (it contradicts some statements in the litterature about the "failure" of 
the post-Newtonian approximation in the regime of the ICO), but we shall present arguments supporting 
it. In particular the 3PN order is probably able to locate the ICO to within an accuracy of the order of 
1% or better (for binary systems with comparable masses). 

2 Post-Newtonian iteration for isolated systems 

Before embarking into strong statements concerning the "accuracy" of the post-Newtonian description of 
binary systems, it is wise to examine in a somewhat general way the well-definiteness of the approximation 
itself. This is especially important in view of the fact that our old lady the post-Newtonian approximation 
has been unjustly accused of being plagued with some apparently inherent difficulties, and furthermore 
which crop up around the 3PN order we are interested in. The problems arise in the general case of regular 
(singularity-free) matter sources. Up to the 2.5PN order the approximation can be worked out without 
problems, and often at the 3PN order the problems can be solved specifically for each case at hands. 
However, it must be admitted that these difficulties, even if appearing at higher approximations, have 
cast doubt in the past on the actual soundness, on the theoretical point of view, of the post- Newtonian 
expansion. They pose the practical question of the reliability of the approximation when comparing the 
theory's predictions with very precise experimental results. Here we discuss the nature of the problems - 
are they purely technical or linked with some fundamental drawback of the approximation? - and outline 
their resolution recently proposed in Ref. [2] ■ 

The first problem is that in higher approximations some divergent Poisson-type integrals appear. 
Recall that the post-Newtonian expansion replaces the resolution of an hyperbolic-like d'Alembertian 
equation by a perturbatively equivalent hierarchy of elliptic-like Poisson equations. Rapidly it is found 
during the post-Newtonian iteration that the right-hand-side of the Poisson equations acquires a non- 
compact support (it is distributed over all space), and that the standard Poisson integral diverges because 
of the bound of the integral at spatial infinity, i.e. r = |x| ^ -l-oo, with t — const. For instance some of 
the potentials occuring at the 2PN order in Chandrasekhar's work [3] are divergent, so the corresponding 
metric is formally infinite. In fact, Kerlick ^ U] showed that the post-Newtonian computation d la 
Chandrasekhar, following the iteration scheme of Anderson and DeCanio 0, can be made well-defined 
up to the 2.5PN order, but that this does not solve the problem at the next 3PN order, which has been 
found to involve some inexorably divergent Poisson integrals. 

These divergencies come from the fact that the post-Newtonian expansion is actually a singular 
perturbation, in the sense that the coefficients of the successive powers of 1/c are not uniformly valid 



2 



in space, since they typically blow up at spatial infinity like some positive powers of r. For instance, 
Rendall T has shown that the post-Newtonian expansion cannot be "asymptotically flat" starting at the 
2PN or 3PN level, depending on the adopted coordinate system. The result is that the Poisson integrals 
are in general badly-behaving at infinity. Trying to solve the post-Newtonian equations by means of the 
standard Poisson integral does not a priori make sense. This does not mean that there are no solution to 
the problem, but simply that the Poisson integral does not constitute the correct solution of the Poisson 
equation in the context of post-Newtonian expansions. So the difficulty is purely of a technical nature, 
and will be solved once we succeed in finding the appropriate solution to the Poisson equation"^. 

To cure the problem of divergencies we have introduced 2 , at any post-Newtonian order, a generalized 
solution of the Poisson equation with non-compact support source, in the form of an appropriate finite 
■part of the usual Poisson integral: namely we regularize the bound at infinity of the Poisson integral by 
means of a process of analytic continuation, analogous to the one already used to regularize the retarded 
integrals in Refs. [51 ll2[ IUlll(J | . At any post-Newtonian order n we have to solve a Poisson equation with 
non-compact-support "source-term" r^"^, where the overbar indicates that we are considering a (formal) 
post-Newtonian expansion, and pick up the coefficient of 1/c" in that expansion (/^, are space-time 
indices). We multiply the source-term by a regularization factor r-^, where r = |x| and B E C, and then 
apply the standard Poisson integral. The result, 

A-'[.-^r](x,0 = -3^X,^ly|».r(y.O. (2) 

is a well-defined integral (i.e. convergent at infinity, r -\-oo) in some region of the i3-complex plane, 
given by 3?(i?) < — a,nax ^ 2, where Omax denotes the maximal power of r in the behaviour of the source- 
term t!^"^ when r — s- -|-oo, i.e. it corresponds to the maximal growth of the source-term at infinity (notice 
that Cmax gets larger and larger when we increase the post-Newtonian order n). Next, we can prove that 
the latter function of B generates a (unique) analytic continuation down to a neighbourhood of the origin 
-6 = 0, except at i? = itself, around which value it admits a Laurent expansion with multiple poles 
up to some finite order. Then, we consider the Laurent expansion of that function when B —^ and 
compute the finite part {J-V), or coefficient of the zero-th power of B, of that expansion. This defines 
our generalized Poisson integral: 

A^i[Tr] =^PA-M^''^r] ■ (3) 

The resulting "generalized" Poisson integral constitutes an appropriate solution of the post- Newtonian 
equation. However this is only a particular solution of the Poisson equation, and the most general 
solution will be the sum of that particular solution and the most general solution of the corresponding 
homogeneous equation. At this stage, considering the post-Newtonian iteration scheme alone, we cannot 
do more and therefore we leave the homogeneous solution unspecified (we can see a posteriori that it is 
associated with radiation-reaction effects). 

The second problem has to do with the near-zone limitation of the post-Newtonian approximation. 
Indeed the post-Newtonian expansion assumes that all retardations r/c are small, so it can be viewed as 
a formal near-zone expansion when r — s- 0, which is valid only in the region surrounding the source that 
is of small extent with respect to the typical wavelength of the emitted radiation: r ^ A (if we locate 
the origin of the coordinates r = inside the source). Therefore, the fact that the coefficients of the 
post-Newtonian expansion blow up at spatial infinity, when r — > -|-c», has nothing to do with the actual 
behaviour of the field at infinity. The serious consequence is that it is not possible, a priori, to implement 
within the post-Newtonian iteration the physical information that the matter system is isolated from the 
rest of the universe. Most importantly, the no-incoming radiation condition, imposed at past null infinity, 
cannot be taken into account, a priori, into the scheme. In a sense the post-Newtonian approximation 

^ The problem is somewhat similar to what happens in Newtonian cosmology. Here we have to solve the Poisson equation 
AU = —iirGp, where the density p of the cosmological fluid is constant all over space: p = p{t). Clearly the Poisson integral 
of a constant density does not make sense, as it diverges at the bound at infinity like the integral J rdr. This nonsensical 
result has occasionally been referred to as the "paradox of Seeliger" . However the problem is solved once we realize that the 
Poisson integral does not constitute the appropriate solution of the Poisson equation in the context of Newtonian cosmology. 
A well-defined solution is simply given by f7 = —^nGpr^. 
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is not "self-supporting" , because it necessitates some information taken from outside its own domain of 
validity. 

The solution of the problem of the near-zone limitation of the post- Newtonian expansion resides in the 
matching of the near-zone field to the exterior field, a solution of the vacuum equations outside the source 
which has been developed in previous work |S] using some post- Minkowskian and multipolar expansions. 
In the case of post-Newtonian sources, the near zone, i.e. r <C A, covers entirely the source, because the 
source's radius itself is such that a <ti X. Thus the near zone overlaps with the exterior zone where the 
multipole expansion is valid. What we do is to impose a matching condition resulting from the numerical 
equality between the multipolar and post-Newtonian fields, verified in the external part of the near-zone 
(say a < r <C A, where a is the size of the source). The matching equation reads 

'M(M^=M(h'"") , (4) 

where A4 symbolizes the multipolar (actually multipolar-post-Minkowskian) series, and as before the 
overbar refers to the post-Newtonian or near-zone expansion. We emphasize that the matching equation 
is physically justified only for post-Newtonian sources, for which the exterior near-zone exists. 

The requirement of matching to the post-Newtonian solution has been shown (in a previous work 
to entirely determine, up to any post-Newtonian order, the multipole moments parametrizing the 
exterior field. In our more recent paper [2| we have proved that the still undetermined homogeneous 
solutions in the post-Newtonian iteration alluded to before are also fully determined by the matching 
equation. These homogeneous solutions are associated with radiation-reaction effects — for instance they 
incorporate the dominant radiation-reaction force at the 2.5PN order, as well as the tail contribution in 
the radiation-reaction force which is known to arise at the 4PN order I12L [T^ l. In conclusion, the 
post-Newtonian expansion of the field inside an isolated system is fully determined (it can be indefinitely 
reiterated, up to any post-Newtonian order) by the matching to the exterior solution satisfying the correct 
boundary condition at infinity — notably the absence of incoming radiation from past null infinity. 

3 Dynamics and radiation of compact binaries 

By equations of motion we mean the explicit expression of the accelerations of the compact bodies in 
terms of the positions and velocities. In Newtonian gravity, writing the equations of motion for a system 
of N particles is trivial; in general relativity, even writing the equations in the case = 2 is difficult. The 
first relativistic term, at the IPN order, was derived by Einstein, Infeld and Hoffmann J3] by means of 
their famous "surface-integral" method, by which the equations of motion are deduced from the vacuum 
field equations, and so which is applicable to any compact objects (they be neutron stars, black holes, 
or, maybe, naked singularities). 

Concerning the 2PN and 2.5PN approximations the result for the equations of binary motion in 
harmonic coordinates was obtained by Damour and Deruelle and collaborators |18l 1151 1161 [T7j . The 
corresponding result for the ADM-Hamiltonian of two particles at the 2PN order was given in Ref. ^H] ■ 
By ADM-Hamiltonian we refer to the Fokker-type Hamiltonian, which is obtained from the matter- 
plus-field Arnowitt-Deser-Misner Hamiltonian by eliminating the field degrees of freedom. Kopeikin 
derived the 2.5PN equations of motion for two extended compact objects. The 2.5PN-accurate harmonic- 
coordinate equations as well as the complete gravitational field (namely the metric g^^i,) were computed 
by Blanchet, Faye and Ponsot and by Pati and Will [221 ■ 

It is important to realize that the 2.5PN equations of motion are known to hold in the case of binary 
systems of strongly self-gravitating bodies. This is via an "effacing" principle for the internal structure 
of the bodies. As a result, the equations depend only on the "Schwarzschild" masses, mi and m2, of the 
compact objects, as has been explicitly verified up to the 2.5PN order |20[ I23|. The 2.5PN equations 
of motion have also been established by Itoh, Futamase and Asada j^|2Sli who use a variant of the 
surface-integral approach valid for compact bodies. 

The present state of the art is the 3PN order. The equations of motion at this order have been worked 
out independently by two groups, by means of different methods, and with equivalent results. On one 
hand, Jaranowski and Schafer |26[ I27|. and Damour, Jaranowski and Schafer |28l I29L ISn| . employ the 
ADM-Hamiltonian formalism of general relativity; on the other hand, Blanchet and Faye |31ll84ir^lH3| . 
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de Andrade, Blanchet and Faye and Blanchet and Iyer [3^1, founding their approach on the post- 
Newtonian iteration initiated in Ref. '21' , compute directly the equations of motion from which they infer 
the Lagrangian (instead of the Hamiltonian) in harmonic coordinates. The end results have been shown 
|29[ 05] to be physically equivalent in the sense that there exists a unique "contact" transformation 
of the dynamical variables, that changes the harmonic-coordinates Lagrangian obtained in Ref. |35| 
into a new Lagrangian, whose Legendre transform coincides exactly with the Hamiltonian given in Ref. 
|28j . The 3PN equations of motion, however, depend on one unspecified numerical coefficient, ^static 
in the ADM-Hamiltonian formalism and A in the harmonic-coordinates approach, which is due to some 
incompleteness of the Hadamard self- field regularization method. This coefficient has been fixed by means 
of a dimensional regularization in Ref. PUj (see Section O). 

So far the status of the post-Newtonian equations of motion is quite satisfying. There is mutual agree- 
ment between all the results obtained by means of different approaches and techniques, whenever it is 
possible to compare them: point-particles described by Dirac delta- functions, extended post-Newtonian 
fluids, surface-integrals methods, mixed post-Minkowskian and post-Newtonian expansions, direct post- 
Newtonian iteration and matching, harmonic coordinates versus ADM-type coordinates, different pro- 
cesses or variants of the regularization of the self field of point-particles. 

Let us remark that the 3PN equations of motion are merely "Newtonian" as regards the radiative 
aspects of the problem, because with that precision they contain the radiation reaction force at only the 
lowest 2.5PN order. From the conservative part of the 3PN dynamics (neglecting the 2.5PN radiation 
reaction) we shall obtain the binary's center-of-mass energy E at the 3PN order. We now want to compute 
the variation of E because of the emission of gravitational radiation. For this purpose we replace the 
knowledge of the radiation reaction force in the local equations of motion by the computation of the total 
energy flux in gravitational waves, say £, and we apply the energy balance equation 

Therefore in our approach the computation of the center-of-mass energy E constitutes only one "half" 
of the solution of the problem. The second "half, that of the computation of the energy flux £, is to be 
carried out by application of a wave-generation formalism. 

Following earliest computations at the IPN level [371 138j (at a time when the post-Newtonian correc- 
tions in C had a purely academic interest), the energy flux of inspiralling compact binaries was completed 
to the 2PN order by Blanchet, Damour and Iyer [21], using the 2PN wave-generation formalism of Ref. 
[H], and, independently, by Will and Wiseman 001, using their own formahsm (see Refs. ^202] for 
joint reports of these calculations). The preceding approximation, 1.5PN, which represents in fact the 
dominant contribution of tails in the wave zone, had been obtained in Refs. 021 by application of a 
formula for tail integrals given in Ref. 0^1 . Higher-order tail effects at the 2.5PN and 3.5PN orders, as 
well as a crucial contribution of tails generated by the tails themselves (the so-called "tails of tails" ) at 
the 3PN order, were obtained by Blanchet gHlllZl- However, unlike the 1.5PN, 2.5PN and 3.5PN orders 
that are entirely composed of tail terms, the 3PN approximation involves also, besides the tails of tails, 
many non-tail contributions coming from the relativistic corrections in the (source) multipole moments 
of the binary. These have been completed by Blanchet, Iyer and Joguet ^HIEHIj based on the expressions 
for the multipole moments given in Ref. |10| . except for one single unknown numerical coefficient, due 
to the use of the Hadamard regularization, which is a combination of the parameter A in the equations 
of motion, and a new parameter coming from the computation of the 3PN quadrupole moment. 

The post-Newtonian flux C is in complete agreement, up to the 3.5PN order, with the result given by 
the very different technique of linear black-hole perturbations, valid in the "test-mass" limit where the 
mass of one of the bodies tends to zero (limit — > 0, where z/ = /x/m). Linear black-hole perturbations, 
triggered by the geodesic motion of a small mass around the black hole, have been applied to this problem 
|50[ I51| . This technique has culminated with the beautiful analytical methods of Sasaki, Tagoshi and 
Tanaka |52l I53[l^ . who solved the problem up to the extremely high 5.5PN order. 



4 Post-Newtonian templates for compact binary inspiral 

The orbital phase of the binary — i.e. the integral of the frequency w(t): 
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m = J Lu{t)dt , (6) 

constitutes the crucial quantity to be monitored (and therefore to be predicted) in the detectors. The 
templates built from the theoretical prediction for the orbital phase should be accurate enough over most 
of the inspiral phase, within the frequency bandwidth of the detectors, with reduced cumulative phase 
lags, so that the phasing errors are not significant when one attempts to extract the values of the binary's 
parameters (essentially the masses and spins) from the data. 

The relevant model for describing the inspiral phase consists of two point-masses moving under their 
mutual gravitational attraction. As a simplification for the theoretical analysis, the orbit of inspiralling 
binaries can be considered to be circular, apart from the gradual inspiral, with a good approximation. 
The templates are based on the energy-balance equation ©, from which one deduces the orbital phase 
as 

where $c and ujc denote the values at the instant of coalescence. The number of gravitational-wave cycles 
left from the current time till the coalescence instant (we consider only the dominant harmonics at twice 
the orbital frequency) is 

J^^-l^—^. (8) 

TT 

It is clear — because of Eq. — that J\f is of the order of the inverse of radiation-reaction effects, 
hence the formal "post-Newtonian" order is M = O (c"*"^). 

As a matter of fact, J\f will be a large number, approximately equal to 1.6 x 10'* in the case of two 
neutron stars between 10 and 1000 Hz (roughly the frequency bandwidth of the detector VIRGO). Data 
analysts have estimated that, in order not to suffer a too severe reduction of signal-to-noise ratio, one 
should monitor the phase evolution with an accuracy comparable to one gravitational-wave cycle (i.e. 
5Af ^ 1), over the whole detector's bandwidth. From a strict post-Newtonian point of view, we see that 
the consequence, since the "Newtonian" number of cycles is formally O {c^^), is that any post-Newtonian 
correction therein that is larger than the order c^^ is expected to contribute to the phase evolution more 
than what is allowed by the previous estimate. Therefore, one expects that in order to construct accurate 
enough templates it will be necessary to compute the orbital phase up to at least the 2.5PN order at 
a minimum. This back-on-the-envelope estimate has been confirmed by measurement-accuracy analyses 
|56|l51llK7ll58| which showed that in advanced generations of detectors the 3PN approximation (or, even 
better, the 3.5PN one) is required in the case of inspiralling neutron star binaries. 

The first ingredient in the theoretical calculation is the binding energy E (in the center-of-mass frame) , 
defined as being the invariant energy associated with the conservative part of the binary's 3PN dynamics 
(we ignore the radiation reaction effect at the 2.5PN order). The center-of-mass energy E is deduced from 
the Hamiltonian in ADM- type coordinates |2f)l 1271 1281 1291 IHII] , or equivalently from the Lagrangian in 
harmonic-coordinates 1321 1331 1341 1351 EE] ■ Restricting our consideration to circular orbits, the energy 
is a function or a single variable, the radial distance r between the two particles in a given coordinate 
system. In fact it is better to express the energy in terms of the frequency w = ^ of the orbital motion, 
or, rather, in terms of the particular frequency-related parameter 



(9) 



The individual masses of the black holes are denoted by mi and m2, and the total mass is M = mi +m2. 
The interest of using the parameter x, instead of some coordinate distance r, is that the energy function 
E{x) then takes an invariant invariant (the same in different coordinate systems). The result we get 
consists of the Newtonian contribution, proportional to x, followed by post-Newtonian corrections up to 
the 3PN order: 
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This expression involves the useful ratio between reduced and total masses: 
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This ratio is interesting because of its range of variation: < v < j, where v = j in the equal-mass 
case and v in the test-mass limit for one of the bodies. The parameter A denotes a point-mass 
regularization-constant and will be discussed in Section [S] 

The second ingredient in this analysis concerns the gravitational-wave luminosity £, calculated, in 
the post-Newtonian approximation, from a wave-generation formalism valid for extended "fluid" systems 
1^1 137| El , and then specialized to binary systems of point-particles ,39,, ^46., J7, .48^ 49^ . The calculation 
takes properly into account the relativistic corrections linked with the description of the source (multipole 
moments), as well as the non-linear effects in the propagation of the waves from the source to the far 
zone. The crucial input of any post-Newtonian computation of the flux is the mass quadrupole moment 
(because the post-Newtonian precision required for the higher multipole moments is smaller). The 3PN 
quadrupole moment for circular binary orbits in a harmonic coordinate system is of the form** 



lij ^ [ Axij + B — Uij 



(12) 



and has been obtained recently by Blanchet, Iyer and Joguet with result 
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where the post-Newtonian ordering parameter (defined in harmonic coordinates) reads 
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Notice the logarithms entering these formulas at the 3PN order. One type involves a constant length 
scale ro coming from the general wave-generation formalism of Refs. 1471 11(1) . and which corresponds 
to some "infrared" cut-off in the bound at infinity of the integrals defining the multipole moments. The 
constant ro should and will cancel out when considering the complete multipole expansion of the field 
exterior to the source. The other type of logarithm contains a different length scale Tq, which should 
rather be viewed as an "ultra-violet" cut-off associated with the singular behaviour of the metric at the 
location of the point-particles; Tq is related to the two constants r\ and r'^ — one for each particles — 
appearing in the 3PN equations of motion of Refs. [311 1341 1331 I36j by 



In Tn 



mi 



In t'y 



m2 



Inr'r, 



(16) 



*We neglect a 2.5PN term and denote e.g. 
the two bodies and r = [x|. 
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iXj — ^Sijr'^, where Xi is the harmonic-coordinate separation between 
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As we know that r[ and r'2 are gauge constants, i.e. they can be ehminated by a gauge transformation at 
the 3PN order, Tq will necessarily disappear from our physical, gauge-invariant, result at the end. Besides 
the harmless constants tq and Tq, there are three point-mass regularization-constants in Eqs. (|13|l - (|14() : 
^, K and C (see Section ISJ. However, we shall see that the energy flux for circular orbits depends on one 
combination only of these constants: 6* = ^ -I- 2k -I- C. 

Through 3.5PN order, the result is decomposed into an "instantaneous" part, i.e. generated solely 
by the multipole moments of the source, and a tail part. What we call here the tail part is in fact 
a complicated sum of "tails", "tail squares", and "tails of tails", as computed in Ref. 47 . We find, 
respectively, 
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The Newtonian result has been factorized out in front. Here C — 0.577 • • • denotes the Euler constant. As 
we can see, the constant ro drops out from the sum of Eqs. H17|) and (|18|l . However, the gauge constant 
Tq does not seem to disappear at this stage, but that is simply due to our use of the post-Newtonian 
parameter 7 defined by H15|l , and which depends via the equation of motion on the choice of harmonic 
coordinates. After substituting the frequency-related parameter x in place of 7 (with the help of the 
3PN equations of motion), we find that Tq does cancel as well — this nicely illustrates the consistency 
between our two computations, in harmonic-coordinates, of the equation of motion on one hand and the 
multipole moments on the other hand. Finally we obtain 



Z = 



32c- 



1247 35 \ , V2 

V X + Attx^' 
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16285 176419 



176 88 
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19897 
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3024 ^ 
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324' 
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1512 



378 



7/2 



(19) 



The last test (but not the least) is that Eq. H19(l is in perfect agreement, in the test-mass limit — > 0, 
with the result of linear black- hole perturbations obtained by Tagoshi and Sasaki ■ 



5 Problem of regular izat ion ambiguities 

A model of structureless point masses is expected to be sufficient to describe the inspiral phase of compact 
binaries. Thus we want to compute the metric (and its gradient needed in the equations of motion) at 
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the 3PN order for a system of two point-like particles. Applying general expressions for the metric valid 
in the case of continuous (smooth) matter sources, we find that most of the integrals become divergent 
at the location of the particles. Consequently we must supplement the calculation by a prescription for 
how to remove (i.e. to regularize) the infinite part of these integrals. 

The "standard" Hadamard regularization yields some ambiguous results for the computation of cer- 
tain integrals at the 3PN order, as Jaranowski and Schafer 26, 27 noticed in their computation of 
the equations of motion within the ADM-Hamiltonian formulation of general relativity. They showed 
that there are two and only two types of ambiguous terms in the 3PN Hamiltonian, which were then 
parametrized by two unknown numerical coefficients Wstatic and Wkinetic- 

Blanchet and Faye |32[I33| . motivated by the previous result, introduced their "improved" Hadamard 
regularization, based on a theory of pseudo- functions and generalized distributional derivatives. This new 
regularization is mathematically well-defined and free of ambiguities; in particular it yields unique results 
for the computation of any of the integrals occuring in the 3PN equations of motion. Unfortunately, 
this regularization turned out to be in a sense incomplete, because it was found |31[ that the 3PN 
equations of motion involve one and only one unknown numerical constant, called A, which cannot be 
determined within the method. The comparison of this result with the work of Jaranowski and Schafer 
|2(il I27| . on the basis of the computation of the invariant energy of binaries moving on circular orbits, 
showed that 



41 

^kinetic = ^ 7 (20) 
11 1987 

^static — ^"3" "§40" • ' ^ 

Therefore, the ambiguity Wkinetic is fixed, while A is equivalent to the other ambiguity Wstatic- The 
value of Wkinctic given by Eq. (|2()|l was recovered by Damour, Jaranowski and Schafer |28', who proved 
that this value is the unique one for which the global Poincare invariance of their formalism is verified. 
By contrast, the harmonic-coordinate conditions preserve the Poincare invariance, and therefore the 
associated equations of motion should be Lorentz-invariant, as was indeed found to be the case by 
Blanchet and Faye |31l I34| , thanks in particular to their use of a Lorentz-invariant regularization |33j 
(hence their determination of Wkinctic)- 

The appearance of one and only one physical undeterminacy A <4> ^static in the equations of motion 
constitutes a quite striking fact, specifically related to the use of an Hadamard- type regularization. 
Mathematically speaking, the presence of A is (probably) related to the fact that it is impossible to 
construct a distributional derivative operator satisfying the Leibniz rule for the derivation of the product. 
The Einstein field equations can be written into many different forms, by operating some terms by parts 
with the help of the Leibniz rule. All these forms are equivalent in the case of regular sources, but they 
become inequivalent for point particles if the derivative operator violates the Leibniz rule. 

In Ref. 60, it has been argued that the numerical value of the parameter ^static could be ~ —9, because 
for such a value some different resummation techniques, viz Pade approximants |58| and effective-one- 
body (EOB) method [21], when they are implemented at the 3PN order, give approximately the same 
numerical result for the location of the last stable circular orbit. Even more, it was suggested [SII that 
ti^static might be precisely equal to Ltj*^^^^^^, with 



^static - - — ^ — • (22) 



47 41 

1 

3 64 

(We have i-iJ static ~ —9.34 • • •.) However, the value of ^static in general relativity has been computed by 
Damour, Jaranowski and Schafer j3(Jj by means of a dimensional regularization, instead of an Hadamard- 
type one, within the ADM-Hamiltonian formalism. Their result is 

1987 

a;.,.,,,^0 ^ ^ = -3080- (23) 

As Damour et al |30| argue, clearing up the ambiguity is made possible by the fact that the dimensional 
regularization, contrary to the Hadamard regularization, respects all the basic properties of the algebraic 
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and differential calculus of ordinary functions^. 

Let us comment tliat the use of a self-field regularization in this problem, it be dimensional or based 
on the Hadamard partie finie, signals a somewhat unsatisfactory situation on the physical point of view, 
because we would like to perform, ideally, a complete calculation valid for extended bodies, taking into 
account the details of the internal structure of the bodies (energy density, pressure, internal velocities, 
etc.). By considering the limit where the radii of the objects tend to zero, one should recover the same 
result as obtained by means of the point-mass regularization, and determine the value of the regularization 
parameter A. Because of considerable difficulties arising at the 3PN order this program has not yet been 
achieved. 

Concerning the 3PN radiation field of two point masses — the second half of the problem, besides the 
3PN equations of motion — , Blanchet, Iyer and Joguet IS' used the (standard) Hadamard regularization 
and found necessary to introduce three additional regularization constants ^, k and C,. However the total 
gravitational-wave flux, in the case of circular orbits, depends only on the linear combination 

6( = ^ 2k + C • (24) 

Furthermore, this 6 comes in at the same level as A, so there is in fact only one unknown constant in the 
flux [see Eqs. (|17() - (|19|l ]. Notice that the improved version of the Hadamard regularization proposed in 
Refs. [H3 OSI should be able, in principle, to fix the value of the constant C,. 

6 Accuracy of the post-Newtonian approximation 

In this section we discuss (some aspects of) the accuracy of the post-Newtonian approximation as re- 
gards the determination of the binary's innermost circular orbit (ICO). The ICO will be defined by the 
minimum, when it exists, of the energy function E{x) given by Eq. Notice that we do not define 

the ICO as an ISCO, i.e. as a point of dynamical general-relativistic unstability. See Ref. [3^ for a 
discussion of the dynamical stability of circular binary orbits at the 3PN order. 

Let us first confront the prediction of the standard (Taylor-based) post-Newtonian approximation 
at the 3PN order for the ICO — as given by the minimum of Eq. Hl()|l — with a recent numerical 
calculation by Gourgoulhon, Grandclement and Bonazzola |61l I62| . These authors obtained numerically 
the energy E{lo) of binary black holes along evolutionary sequences of equilibrium configurations under 
the assumptions of conformal flatness for the spatial metric and of exactly circular orbits. The latter 
restriction is implemented by requiring the existence of an "helical" Killing vector, time-like inside the 
light cylinder associated with the circular motion and space-like outside. The numerical calculation 

I52| has been performed in the case of corotating black holes, which are spinning with the orbital 
angular velocity oj. For the comparison we must therefore include within the post-Newtonian formalism 
yielding Eq. (|10|l the effects of spins appropriate to two Kerr black holes rotating at the orbital rate uj. 

The Figure n (issued from Ref. fH3l) presents the post-Newtonian results for Eico in the case of 
irrotational and corotational binaries. The points indicated by IPN, 2PN and 3PN are defined from 
the obvious truncation of Eq. The points IPN'^"™*, 2PN™"* and SPN'^"™* take into account the 

spin effects of corotational binaries and are computed in Ref. |63| . Notice that the irrotational and 
corotational configurations differ only from the 2PN order. As we can see the 3PN points, and even the 
2PN ones, are rather close to the numerical value (marked by an asterisk in Figure As expected, the 
best agreement is for the 3PN approximation and in the case of corotation: i.e. the point SPN'^"™*. So 
our first conclusion is that the location of the ICO as computed by numerical relativity, under the helical- 
symmetry approximation is in good agreement with post-Newtonian predictions. This constitutes 
an appreciable improvement of the previous situation, because we recall that the earlier estimates of the 
ICO in post-Newtonian theory ^64, and numerical relativity |65l I66| strongly disagree with each other, 
and do not match with the present 3PN results (see Ref. [221 for further discussion). 

Our second conclusion comes from the fact that the 2PN and 3PN values in Figure — either for 
irrotational or corotational binaries — are so close to each other. Indeed it seems that the 3PN points 

^Note that the result 1231 is very different from lu*^^^^'^^ given by Eq. I22i : this suggests that the resummation techniques 
(Pade approximants and EOB method), though they are designed to "accelerate" the convergence of the post-Newtonian 
series, do not in fact converge toward the same exact solution (or, at least, not as fast as expected). See Section 171 for 
discussion on this point. 
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Figure 1: The center-of-mass energy Eico versus wico in the equal- mass case {f = j). The asterisk 
marks the resuh calculated by numerical relativity |^|S2]- The points indicated by IPN, 2PN and 3PN 
correspond to irrotational binaries, while the points denoted by IPN'^™"', 2PN'^°™* and SPN'^™"' describe 
corotational binaries. Both 3PN are 3PN'^°''°* are shown for Wstatic = 0. 

are useful only to confirm the result given by the 2PN ones. This is a quite satisfying state of affairs 
because it indicates that the post-Newtonian approximation converges very well (in the sense of Cauchy's 
criterion). Figure ^ shows that the post-Newtonian approximation is well qualified to accurately locate 
the ICO (of course for this purpose one must go to the 3PN order — the IPN approximation is clearly 
not accurate enough). 

Let us elaborate more on this point. First we make a few order-of-magnitude estimates. At the 
location of the ICO we find (see Figure that the frequency-related parameter x defined by Eq. © 
is approximately of the order of 20%. Therefore, we might a priori expect that the contribution of the 
IPN approximation to the energy at the point of the ICO should be of the order of 20%. For the present 
discussion we take the pessimistic view that the order of magnitude of an approximation represents also 
the order of magnitude of the higher-order terms that are neglected. We see that the IPN approximation 
should yield a rather poor estimate of the "exact" result, but this is quite normal at this very relativistic 
point where the orbital velocity is - ~ y/x ~ 50%. By the same argument we infer that the 2PN 
approximation should do much better, with fractional errors of the order of ~ 5%, while the 3PN 
approximation will be even better, with the precision x^ ~ 1%. 

The simple order-of-magnitude estimate suggests therefore that the 3PN order should be close to the 
"exact" solution for the ICO to within 1% of fractional accuracy. We think that this is very good, and 
we should even remember that this estimate is pessimistic, because we could reasonably expect that the 
neglected higher-order approximations, 4PN and so on, are in fact much smaller numerically (e.g. of the 
order of x"^ ~ 0.2%). But let us keep for the present discussion the 1% guess for the accuracy of the 3PN 
approximation. 

Now the previous estimates make sense only if the numerical values of the post-Newtonian coefficients 
in Eq. H10|l are roughly of the order of one. If this is not the case, and if the coefficients increase 
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dangerously with the post-Newtonian order n, one sees that the post-Newtonian approximation might in 
fact be very bad. So let us look at the values of the coefficients of the IPN, 2PN and 3PN approximations 
in Eq. (QUI), say 



ai{v) = 
03(1^) = 



3 jy_ 

i ^ 12 ' 
27 19 

675 



V 

"""24' 
209323 
4032 " 



205 2 110, 



155 
"96"' 



35 
5184^ 



(25) 
(26) 
(27) 



We present in Tabled the values of these coefficients in the test-mass limit v 
case 1/ = i when the ambiguity parameter takes the "uncorrect" value Ci;*^, 
correct one i^-'static = predicted by general relativity j3U| . 



0, and in the equal-mass 
^ [see Eq. and the 





Newtonian 


a\{v) 


a2(y) 




y = 


1 


-0.75 


-3.37 


-10.55 


1 


Static - -9-34 1 


-0.77 


-2.78 


-8.75 



0.77 



-2.78 



-0.97 



Table 1: Sequence of coefficients of the post-Newtonian series composing the energy fmiction H25|l - (|27|l . 

Our first comment is that when v — ^ there is an increase of the coefficients by roughly a factor 3 at 
each step. This behaviour is fairly easy to understand. It comes from the existence in the Schwarzshild 
metric of the famous light-ring orbit: a geodesies of photon which is a circular orbit located at i? = 3M in 
Schwarzshild coordinates. As a result, the energy of the test particle in the Schwarzshild metric exhibits 
a singularity at xught-ring = 5, 

The consequence is that the radius of convergence of the post-Newtonian series is ^ , so the post-Newtonian 
coefficients increase by a factor ^ 3. So the post-Newtonian series is not very accurate in the case v = 
This fact has motivated several statements in the literature (see e.g. ESI EZl EH ) j according to which 
the post-Newtonian approximation would be "poorly convergent" , or that there should be a "fundamental 
breakdown" of its validity in the regime of the ICO. This is indeed true in the test-mass limit (v — *■ 0), 
where the post- Newtonian series converges slowly^. 

On the other hand, what happens in the equal- mass case {v = j)? When v — j and we have the 
value ^static notice that the coefficients increase approximately like in the test- mass case v = {). This 
indicates that the gravitational interaction in the case of Wgtatic looks like that in a one-body problem. 
We shall say that when the post-Newtonian coefficients rapidly increase with the order of approximation, 
the interaction is "Schwarzschild-like" . 

Now when v — ^ and the ambiguity parameter takes the correct value (jJstatic = 0, we see that the 
3PN coefficient a3(|) is of the order of minus one instead of ^ —10. We think that this strongly suggests, 
unless 3PN happens to be quite accidental, that the post-Newtonian coefficients in general relativity do 
not increase very much with n, and stay rather of the order of one. This is interesting as it indicates that 
the actual general-relativistic two-body interaction is not Schwarzschild-like. 

It is impossible of course to be very confident about the validity of the previous statement because 
we know only the coefficients up to the 3PN order. Any tentative conclusion based on the 3PN order 

®Let us remark that this negative conclusion does not matter: indeed we shall never use the post-Newtonian approxi- 
mation in the case f — > simply because we know the exact result which is given by Eq. 1101 . The exact result for the 
radiation field is known as well, albeit numerically only 1511 1571 . Therefore we should not worry too much about the poor 
convergence of the post-Newtonian series in the test-mass limit. The post-Newtonian method is useless and even one might 
say irrelevant when considering the motion of a test particle around a Schwarzschild black hole. 
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can be "falsified" when we obtain the next 4PN order. Nevertheless, we feel that the mere fact that 
a3(i) = —0.97 in Table^is sufficient to motivate our (tentative) conclusion that the field of two bodies is 
more complicated than the Schwarzschild space-time. This conclusion is in accordance with the present 
author's respectfulness of the complexity of the Einstein field equations. 

The nice consequence is that because the post-Newtonian coefficients when 1^ = 2 stay of the order 
of one, the standard post-Newtonian approach, based on the standard Taylor approximants, is probably 
very accurate. The post-Newtonian series seems to "converge well" , with a "convergence radius" of the 
order of one^. A convincing support of this view is provided by Figure ^ itself . The order-of-magnitude 
estimates we did at the beginning of this Section are probably correct. In particular the 3PN order should 
be close to the "exact" solution even in the regime of the ICO. 

It is also interesting to look at the numerical values of the post-Newtonian coefficients in the total 
flux £ obtained in Eq. H19|l . For this case we do not have (like for the ICO) a clear point at which we 
can compare with a result of numerical relativity. Furthermore things concerning the convergence of the 
post-Newtonian series are less clear than with the energy function. A possible reason for that could be 
that the flux, in contrast to the energy, includes the contribution of tails at the 1.5PN, 2.5PN and 3.5PN 
orders, and even of tails-of-tails at the 3PN order. 

Newtonian 63/2(1^) feC^^) br^/2M bsiu) 

u = Q 1 1237 '^im -38.29 128.85 -101.51 

u = \ A = -ijjl^ e = i -4.44 12.57 -0.10 -55.80 115.26 047 

Table 2: Post-Newtonian coefficients in the flux function (|19l) . The coefficient 63(1^) contains a log-term 
that we have evaluated at xico ~ 0.2. [The value 6* = is taken to be indicative.] 



7 On Schwarzschild-like templates for binary inspiral 

Let us finally comment about a possible implication of our conclusion as regards the validity of the 
so-called post-Newtonian resummation techniques, i.e. Pade approximants (58t i60) . which aim at "accel- 
erating" the convergence of the post- Newtonian series in the pre-coalescence stage, and effective-one-body 
(EOB) methods [nnHSIIli which attempt at describing the late stage of the coalescence of two black holes. 
These techniques are based on the idea that the gravitational two-body interaction is a "deformation" 
— with V < ^ being the deformation parameter — of the Schwarzschild space-time. 

The Pade approximants are valuable tools for giving accurate representations of functions having some 
singularities. In the problem at hands they would be justified if the "exact" expression of the energy 
[whose 3PN expansion is given by Eqs. (|25|l - 127|l ] would admit some Schwarzschild-like features, and in 
particular a light-ring singularity. In the Schwarzschild case the Pade series converges rapidly toward the 
solution [SS] : the Pade constructed only from the 2PN approximation of the energy — keeping only ai (0) 
and 02(0) in Eqs. ^ — already coincide with the exact result given by Eq. (jSJ. On the other hand, the 
EOB method maps the post-Newtonian two-body dynamics (at the 2PN or 3PN order) on the geodesic 
motion on some effective metric which happens to be a j/-deformation of the Schwarzschild space-time. 
In the EOB method the effective metric looks like Schwarzschild by definition, and we might expect the 
two-body interaction to own some Schwarzschild-like features. 

Our comment is that the validity of these resummation techniques (and the templates for binary 
inspiral built on them) is questionable, because as we have seen in Section|Blthe value ^static = suggests 
that most probably the two-body interaction is not Schwarzschild-like. In particular, there does not seem 
to exist something like a light-ring orbit which would be a deformation of the Schwarzschild one. 

To test the previous idea, let us come back to the Schwarzschild limit, where we observed that the 
radius of convergence of the post-Newtonian series is given by the light-ring singularity at the value i. 
Now the radius of convergence is given by d'Alembert's criterion as the limit when n — )■ +00 of the ratio 

'^Actually, the post-Newtonian series could be only asymptotic (hence divergent), but nevertheless it should give good 
results provided that the series is truncated near some optimal order of approximation. In this discussion we assume that 
the 3PN order is not too far from that optimum. 
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of successive coefficients, i.e. an-i/cin- We might therefore compute the hght-ring orbit by investigating 
the hmit 

^tcst Y 

a^iight-ring = lim = q ■ (29) 

To test for the possible existence of a hght-ring singularity in the case of comparable masses, we consider 
the ratio between the two highest known post-Newtonian coefficients, that are a2{y) and a^{v). This 
ratio will give an estimate of the "light-ring" singularity for non-infinitesimal mass ratios, 

a;iight-ring -r-r ■ 30) 

Using the values given in Table ^ we obtain for equal masses in the case of the "wrong" ambiguity 
parameter w*^^^;^ ~ -9.34, 

2:iight-ring(|:, ^static) ^ 0-32 . (31) 

As we see there seems to be in this case a (pseudo-)light-ring orbit which is a small deformation of the 
Schwarzschild light-ring orbit given by Eq. (|29|l . In this Schwarzschild-like situation, we expect that 
Pade approximants and EOB-type methods to be appropriate. 

But, in the case of general relativity (wstatic = 0), we obtain a drastically different result: 

a^iight-ring(iG.i?.)- 2.86. (32) 

If we believe the correctness of this estimate we must conclude that there is in fact no notion of a light-ring 
orbit in the real two-body case. Or, one might say (pictorially speaking) that the light-ring orbit gets 
hidden inside the horizon of the final black-hole formed by coalescence. If we consider the ratio between 
the IPN and 2PN coefficients (instead of between 2PN and 3PN), we get the value ~ 0.28 instead of 
Eq. (|32|l . So at the 2PN order the field seems to admit a light ring, while at the 3PN order it does 
not. This reinforces our idea that it is meaningless (with our present 3PN-based knowledge, and untill 
fuller information is available), to assume the existence of a light-ring singularity in the equal-mass case. 
Our expectation, therefore, is that the conditions under which the Pade and EOB methods should be 
legitimate might in fact not be fulfilled. 

This doubt is confirmed by the finding of Ref. |]60j (already alluded to above) that in the case of the 
"wrong" ambiguity parameter w*tatic — —9-34 the Pade approximants and the EOB method at the 3PN 
order give the same result for the ICO. From the previous discussion we see that this agreement is to be 
expected because a deformed (pseudo-)light-ring singularity seems to exist with w*tatic. By contrast, in 
the case of general relativity, where Wstatic = 0, the Pade and EOB methods give quite different results (c/ 
the figure 2 in Ref. |60| ) . Another confirmation comes from the light-ring singularity which is determined 
from the Pade approximants at the 2PN order (see Eq. (3.22) in Ref. |SH1) as a;[[ght™„g ( ^ ) ~ 0.44. 
This value is rather close to Eq. H31|) but strongly disagrees with Eq. H32|l . Our explanation is that the 
Pade series might converge toward a theory having ^static = 'jJltsxtic ^^'i therefore which is different from 
general relativity. 
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